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ABSTRACT 

We investigate whether nuclear star clusters and supermassive black holes follow a common set of 
mass scaling relations with their host galaxy's properties, and hence can be considered to form a single 
class of central massive object. We have compiled a large sample of galaxies with measured nuclear 
star cluster masses and host galaxy properties from the literature and fit log-linear scaling relations. 
We find that nuclear star cluster mass, Mnc, correlates most tightly with the host galaxy's velocity 
dispersion: log Mnc = (2.11 ± 0.31) log(<r/54) + (6.63 ± 0.09), but has a slope dramatically shallower 
than the relation defined by supermassive black holes. We find that the nuclear star cluster mass 
relations involving host galaxy (and spheroid) luminosity and stellar and dynamical mass, intercept 
with but are in general shallower than the corresponding black hole scaling relations. In particular 
Mnc Mq^ 5 ;^ 15 ; the nuclear cluster mass is not a constant fraction of its host galaxy or spheroid 
mass. We conclude that nuclear stellar clusters and supermassive black holes do not form a single 
family of central massive objects. 

Subject headings: galaxies: dwarf — galaxies: fundamental parameters galaxies: kinematics and 
dynamics — galaxies: nuclei — galaxies: star clusters — galaxies: structure — 



1. INTRODUCTION 

Central massive objects (CMOs) are a common feature 
in galaxies across the Hubble sequence. CMOs take 
the form of either a supermassive black hole (SMBH) 
or a compact stellar structure such as a nuclear stellar 
cluster (NC) or nuclear stellar disk (ND). The masses 
of SMBHs have been shown to correlate with a range 
of host gala xy properties including: stellar velocity 
dispe rsi on, a (iFerrarese fc MerrittlfeOOOt iGebhardt et al.l 



200 0t iGraham et al 
([Graham et alj 120011 



2011 ), stellar c oncen tration 
IGraham fc Driven [20071) ; dy- 
namical mass, Md V n oc a R (iMagorrian et all 119981: 
Marconi fc Huntl 120031: lHaring fc Rbd 120041: IGraham 



2012a) 



2012aD : and luminosity, Lg D h ([Kormendv fc Richstone 



Marconi fc Hun 



t 2003). 



Following the discovery that the luminosity of stellar 
CMOs co rrelates with t hat o f their host bulge in disk 
galaxies (|Balcells et a ll [200l [20071 hereaft er BGP071 
and elliptical galaxies ( Graha m fc Guzman! 120031) . the 
masses of stellar CMOs have also been shown to correlate 
with their host galaxy properties. NC mass has been re- 
ported to correlate with, for early-type galaxies, the host 
galaxy's luminosity, Lp a i and dynamical mass , as given 
by Mcai.dyn oc cr 2 -R 0i Gai (jFerrarese et al.|[2~006al hereafter 
F06). Related correlations have also b een reported with 
the host spheroid's: luminosity, L Sph (jWehner fc Harris! 
[2001 hereafter WH06); stellar mass, M Sph ,* (BGP07); 
dynamical mass, M pph.Hvn (WH06, BGP07) and velocity 
dispersion, a (F06, Graham; 120123). 

These scaling relations are physically interesting be- 
cause they relate objects on very different scales: the 
gravitational sphere of influence of a SMBH is typically 
less than 0.1 per cent of its host galaxy's effective radius, 
R e . This connection is thought to be driven by feed- 
back processes from the CMO (e.g ISilk fc Reesl 119981: 
iCroton et all 120061 : iBooth fc Schavell2009D . but may in- 



stead be related by t he initial central stella r density 
of the host spheroid (| Graham fc Driver! I2007T ) . While 
most studies have focused on feedback from black holes, 
analogous mechanisms driven by nuclear stellar clus- 
ters have been hypothesized ([McLaug hlin et al.1 120061 : 
iMcQuillin fc McLaughlinll2012fi . One potential problem 
with these momentum-conserving feedback arguments, 
as constructed, is that they predict a slope of 4 for both 
the Mbh-c and Mnc-c relations, whereas the obser- 
vations now suggest a sl ope of 5 (jFerrarese fc Merrittl 
l2000t [Graham et al.ll201ll ) and somewhere between 1 and 
2 ([Graham! I2012bf ). respectively. It should however be 
noted that the a term in the models relates to that of 
the dark matter halo rather than the stars, and as such 
they may not be appropriate for comparison. 

F06 and WH06 have argued that SMBHs and NCs fol- 
low a single common scaling relation with Md yn (though 
not with other host galaxy properties). Other investi- 
gations have reached different conclusions, for example 
BGP07 find that NCs do not fall onto the linear relation 
defined by massive central black holes, and conclude that 
any CMO-bulge mass relation that encompasses both 
central black holes and nuclear star clusters must not be 
log-linear. 

F06 have reported that the Mnc - ' 7 relation has a slope 
which is consistent with the Mbh~c relation. How- 
ever expa nding upon th e Mqm q -Q" diagram from Fig- 
ure 8 of IGraham et all ([201 lL IGraham! (j2012bD has 
reported that Mnc k o' 1 to a 2 , whereas Mrh oc 
q 5 for non-barred g al axie^ (jFerrarese fc Merrittl 2000: 
IGraham et all 120111 ). iLeigh et al.l (|2012D alsoreport a 
significantly flatter slope of Mnc oc cr 2 - 73 for a sample of 

1 Barred galaxies tend to have higher velocity dispersions than 
given by the Mbh -0- r elation defined by non-barred galaxies 



(Graha m 2008; Hu 20Q8|). As such, the classical (i.e. all galaxy 
types) Mbh — cr relation has a slope ~ 6. 
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NCs. 

The situation is complicated still further by a blur- 
ring of the division between galaxies containing SMBHs 
or NCs. Since F06, and WH06 who initially found a 
clear division in mass between galaxies hosting a SMBH 
(with Ms p h,dyn > 5 x 10 9 M Q ) and galaxies hosting a NC 
(with Msph.dyn < 5 x 10 9 M Q ), an increasing number of 
galaxie s that host both a SMBH and a NC have been 
found (fGraham fc Driver] l2007t iGonzalez Delgadp et al.l 
l200l ISeth et all 120080 . IGraham fe Spitlerl (120091 here- 
after GS09) observed a transition region from 10 8 < 
M S ph,*/M Q < 10 10 where both types o f nuclei coex- 
ist (see also iNeumaver fe Walcherl 12012). These find- 
ings raise the question of how the combined CMO mass, 
Mbh + Mnc, may scale with the host galaxy properties, 
t hough a larger sa mple of such objects is desired. 

IGraham! (]2012bD updated the Mcmo _ c diagram, first 
published by F06, using an expanded sample of galax- 
ies with directly measured SMBH masses which also in- 
cluded 13 galaxies with both a NC and SMBH. Here 
we re-examine and update the Mqmo versus (i) velocity 
dispersion, (ii) 5-band galaxy magnitude and (iii) dy- 
namical mass diagrams from F06. In addition to the 
abov e sample expan s ion, w e incorporate the NC data set 
from iBalcells et al.l (|2007l ) . We also construct another 
three diagrams involving Mqmo and: if -band luminos- 
ity; total stellar mass, Mc a i,*; and spheroid stellar mass, 
Msph,*- 

Collectively our data represents the largest sample to 
date of NC and host galaxy properties. We have used this 
to investigate their scaling relations and whether they are 
consistent with those for SMBHs. Our sample and data 
are more fully described in Section [2] while in Section [3] 
we present a range of NC and SMBH scaling relations. 
In Section |4] we discuss whether our results support the 
idea of a single common scaling relation for CMOs and 
the implications of our results on a common formation 
mechanism for SMBHs and NCs. Finally, we present a 
summary of our conclusions at the end of Section [4] 

2. SAMPLE AND DATA 

We constructed our sample of nuclear stellar objects 
by combining the data from F06 (51 objects). BGP07 
(17 objects) and GS09 (16 objects). iGrahami (|2012bfi 
added a further 3 objects to the GS09 sample for a total 
of 19 objects, 15 with both a NC and a SMBH and a 
further 4 objects with a NC and only an upper limit on 
M BH - NGC7457 appears in both the BGP07 and GS09 
samples, however the GS09 nuclear cluster properties are 
taken directly from BGP07. We eliminate this duplicate 
galaxy, reducing our sample by one galaxy. This gives 
a final sample of 86 objects with measured Mcmo- The 
observed and derived properties of the nuclear star clus- 
ters and their host galaxies are described in full in the 
following sections, and are presented in Tabled] 

2.1. Nuclear stellar masses 

GS09 provide stellar masses for their nuclei, whereas 
F06 and BGP07 tabulate only magnitudes. For the F06 
objects we derived nuclear stellar object masses following 
F06. We multiplied the total CMO g-band magnitude by 
a mass-to-light ratio, M/L, determined from th e singl e 
stellar population models of Bruzual fc Char iot (2003), 



using the nuclear cluster colors given in lCote et al.l (|2006l ) 
and a stellar population age t — 5 Gyrs. For the BGP07 
objects we derived masses following BGP07, by multiply- 
ing the total CMO JT- band magnitude by M/L K = 0.8 
(jBell fc de Jongl 120011 ) based on typical colors of the 
bulge population. The uncertainties on the nuclear ob- 
ject masses for the F06, BGP07 and GS09 data are given 
by the respective authors as 45 per cent, 33 per cent and 
a factor of 2 respectively. In passing we note that if 
NCs are related to ultra compact dwarf galaxies (e.g. 
iKroupa et al.l I20I0I ) they ma y have a high stell ar M/L 
due to either a bottom-heav y dMieske et al. 2008) or to 



or top- 
, 2009). 



heavy initial mass function (jDabringhausen et al 

BGP07 distinguish between extended nuclear compo- 
nents (11 objects) and unresolved nuclear components 
(also 11 objects - 5 galaxies contain both a resolved 
and unresolved nuclear component), finding that the ex- 
tended components are well fit with an exponential pro- 
file and are thus likely to be nuclear disks (or possibly 
nuclear bars), whereas the unresolved components are 
probably nuclear star clusters. They revealed that the 
disks and clusters follow quite different relations, in terms 
of, for example, how the nuclear disk luminosity scales 
with host galaxy a. For this reason it is important to 
distinguish between nuclear disks and clusters when ex- 
amining the scaling relations of nuclear objects with their 
host galaxy. 

F06 identified three of their objects as containing 
small-scale stellar disks. Based on their published HST 
surface photometry we identify one further object, NGC 
4550 (VCC 1619) as likely containing a nuclear stellar 
disk. These four nuclear disks are also the most ex- 
tended nuclear components in the F06 sample, having 
half-light radii ranging from 26 to 63 pc (for comparison, 
the mean half-light radii for the F06 nuclear objects is 
4 pc). This provides final samples of 76 and 15 nuclear 
clusters and nuclear disks, respectively, from the galaxy 
samples of F06, BGP07 and GS09. 5 galaxies contain 
both a nuclear star cluster and a nuclear disk, thus 86 
unique galaxies with a stellar CMO. 

2.2. Host galaxy and spheroid properties 

We we were not able to obtain every galaxy and 
spheroid property for every object from the literature. 
The number of objects for which we were able to ob- 
tain a given property is indicated in the final row of 
Table [H Velocity di s persion s were obtained from F06, 
lFalconT3arroso et al l (f2 002. for the BGP07 galaxies) 



and IGraham et al.1 ( 2011L for the GS09 galaxies), giv- 



ing 51/76 nuclear star cluster galaxies and 15/15 nuclear 
disk galaxies having measured a. The velocity disper- 
sions were measured in inhomogeneous apertures and we 
do not attempt to correct the measurements to a com- 
mon aperture here. The F06 a values were measured 
from long-sl it observations within aJ_R e aperture. The 
values from iFalcon-Ba rroso et al.l (|2002T ) were obtained 
within a central 1.1 arcsec 2 aperture and corrected to a 
'standard' aperture defined to be equivalent to a circular 
aperture with radi us 1.7 arcsec at the di stance of Coma 
(as establish ed bv | Jorgensen et al.l fl9 95) . The a values 
presented in IGraham et al.l ( 20111) were originally drawn 
from the HyperLeda database (|Paturel et al.| [2003)Fl and 

2 http://leda.univ-lyonl.fr 
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TABLE 1 

Nuclear star cluster, nuclear stellar disc and host galaxy properties 
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Q Q 

y.y 


1 n n 

1U.U 


Q /17 

y.4 ( 


Q Q 

y.y 



y.y 


Q fi 

y.o 


n in 

y. iu 


100 

1U.U 




100 

1U.U 


Q 77 

y. / 1 


Q fi 

y.o 






1 fi 

1U.O 


y.o 


y.o 


1 Q 
1U.O 


Q fi7 

y.o 1 


Q fi 

y.o 




y.o 


1 M 
1U.01 


Q zl 

y.4 


Q A 

y.4 




in n 

iu.y 


Q Q 

y.o 




Q A 

y .4 


1 ^fi 

lU.OO 


Q A 

y.4 


Q A 

y.4 


Q 7 

y. 1 


1 fii 

lU.Ol 


y.4 


Q A 

y.4 




1 1 Q 

11.0 



y.z 



y.z 


Q Q 

y.o 


1 

iu.yo 


Q Q 

y.o 




Q Q 

y.o 








O A 

y .4 


1 1 n^ 
11. UO 


O Q 

y.o 


O Q 

y.o 


Q 9 

y .z 


1 1 0^ 
1 1 .UO 


Q 

y.z 


Q 

y.z 




1 1 Q 
11.0 


Q "1 

y. 1 


Q 1 

y . 1 


y.o 


11 1 fi 
1 1 . lo 


Q 
y.z 


Q 
y.z 


9 
y .z 


I 1 9/1 

I I .Z4 



y.z 



y.z 


Q Q 

y.o 


1 9 1 Q 
1Z. 10 


fi fi 
0.0 




Q fi 

y.o 


11 

1 1 .OO 


Q 1 

y. 1 


Q 1 
y . 1 


Q fi 

y.o 


I 1 Qfi 

I I .Oo 


Q 1 
y. 1 


Q 1 
y . 1 




1 O 1 7 
1Z. 1 / 


fi 7 
0. * 


ft 7 
0. ( 




10 C1 
1Z.01 


c 7 
0. / 


fi 7 
0. ( 




1 00 
1Z.ZZ 


fi fi 
0.0 


fi fi 
0.0 




1 O QI 
1 Z . 1 


ft 7 
0. / 


Q 7 
O. i 




1 1 Ofi 

11. yes 


ft 
o.y 


fi O 

o.y 


9 
y .z 


1 n /io 
iu.4y 


y.o 


^ 
y.o 


9.3 


12.82 


8.5 


8.5 




11.66 


9.0 


9.0 




1 1 Q 
ll.O 


o n 
y.u 


n 
y.u 




12.92 


8.6 


8.6 




1 fii 
IZ.Ol 


ft ^ 
O.O 


fi ^ 
O.O 




11.83 


8.9 






1 7 A 
1Z. / 4 


ft fi 
O.O 


fi fi 
O.O 




1 Q Ofi 
10. Zo 


ft fi 
O.O 


fi fi 
O.O 




1 zll 
1Z.41 


ft 7 
O. ( 


fi 7 
0. 1 


y.o 









y.y 


1 O Q/l 
1Z.04 


ft 7 
O. / 


fi 7 
0. i 




1 Q Ofi 
10. ZO 


ft zl 
0.4 


fi A 
0.4 




10. OU 


ft Q 
O.O 


fi Q 
O.O 


8.9 


5.09 


9.0 


8.41 


11.0 


6.24 


10.8 


10.5 


11.7 


6.31 


11.3 


11.2 


10.7 


9.51 


10.1 


9.6 


11.0 


5.88 


10.9 


10.9 


10.5 


6.75 


10.6 


10.2 


10.6 


7.58 


10.4 


10.0 




10.0 


8.5 


7.5 


11.1 


6.37 


11.1 


11.1 




4.10 


9.5 


8.2 


8.5 


5.59 


8.9 


8.9 



16 



17.2401 


8.06 


15.2801 


8.72 


17.4701 


7.95 



9801 



8.09 



2.80 
2.87 
1.42 
3.40 
3.28 
2.23 
2.29 
3.11 
2.30 
2.55 
1.10 
2.44 
1.39 
1.52 
1.91 
1.41 
1.99 
1.64 
1.33 
1.01 
1.12 
1.93 
1.27 
2.24 
1.33 
1.12 
1.99 
1.07 
1.31 
1.19 
0.85 
0.92 
1.10 
1.14 
1.61 
1.15 
2.46 
2.07 
1.12 
1.23 
1.33 
1.08 
0.75 
1.08 
1.15 
2.09 
1.26 
1.04 
1.72 
0.96 
1.37 



6801 
3701 



3701 
6101 



72M 
204 W 
329 W 
162 W 
252 W 
148W 
178M 

17l'[k] 
37M 
20M 
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Table 1 continued 



Galaxy 


Type 


m s 


m-M 


m N c 


M N c 
log 


m N D 


M ND 
log 


M/L 




Re 


M d y„ 

log 


m K 


M G al 

log 


M Sp h 
log 






mag 


mag 


mag 


M 


mag 


M 




km s _1 


arcsec 


M Q 


mag 


M 


M 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


N3621 


Sd 


9.93l b l 


29.16l f l 




7.00l f l 










27l b l 




6.60 


10.2 


8.2 


N4041 


Sbc 


11.80W 


31.78M 




7.46M 














8.41 


10.6 


8.8 


V1254 


dE 


15.02W 


31.09W 




7.04M 








3lM 


4W 


9.0 


11.70 


9.5 


9.5 


N1300 


SBbc 


10.98W 


31.58W 




7.94W 








229 W 


33W 




7.56 


10.8 




N2549 


SBO 


11.91W 


30.45W 




7.04W 








144M 


6 [b] 


10.4 


8.05 


10.2 




N3585 


SO 


10.60W 


31. 45^ 




6.60W 








206 M 


12W 


11.2 


6.70 


11.1 




Milky 

Way 

N5326 




-0. fo L J 






i .4o L J 








1 nnfkl 
1UU L J 


v>ooo4 L J 




Q 97 


10 7 


in 1 


Sa 


12.83W 


32.68l h ) 






13.47W 


8.90 


0.80 


164W 


12W 


11.3 


8.88 


10.7 


10.4 


N5389 


SO 


12.79W 


32.09l h ' 






15.02W 


8.04 


0.80 


114W 


6 [b] 


10.6 


8.62 


10.6 


10.1 


N5422 


so 


12.75W 


32.02W 


17.23W 


7.13 






0.80 


160W 






8.76 


10.5 


10.0 


N5443 


SBb 


13.29H 


32.17W 






16.18( h ] 


7.61 


0.80 


76W 






9.04 


10.5 


9.5 


N5475 


Sa 


13.37W 


31.95W 


16.95W 


7.22 


13.25l h l 


8.70 


0.80 


9lW 






9.40 


10.2 


9.6 


N5587 


so 


13.74W 


32.46W 






14.23W 


8.51 


0.80 


93m 






9.68 


10.3 


9.3 


N5689 


SO 


12.65W 


32.41W 






13.13t h ] 


8.93 


0.80 


143W 






8.40 


10.8 


10.1 


N5707 


Sab 


13.23W 


32.46W 


16.30W 


7.68 


12.69W 


9.12 


0.80 


14lW 










10.0 


N5719 


SBa 


13.14M 


31.82W 






15.18W 


7.87 


0.80 


108M 


lQ[b] 


10.7 


8.23 


10.7 


10.4 


N5746 


SBb 


11.12W 


31.80W 


17.08W 


7.10 






0.80 


139W 


25W 




6.88 


11.2 


10.1 


N5838 


E/SO 


11.69W 


31.31W 


15.18W 


7.67 


11.15 W 


9.28 


0.80 


255W 


8 M 


11.2 


7.58 


10.7 


10.3 


N5854 


SO 


12.48 [bl 


31.84W 


14.94W 


7.98 


12.70 [h] 


8.87 


0.80 


97W 


5 [b] 


10.3 


8.82 


10.4 


9.9 


N5879 


Sbc 


12.17W 


30.70W 


17.42W 


6.53 






0.80 


58W 


9 [b] 




8.79 


10.0 


9.0 


N6010 


SO 


12.69W 


32.07W 


16.29W 


7.53 






0.80 


144W 






8.93 


10.5 


9.9 


N6504 


Sa 


13.06W 


33.96W 


17.25W 


7.90 






0.80 


185 W 






9.22 


11.1 


10.7 


N7457 


SO 


11.87( b ] 


30.68W 


14.97W 


7.50 


14.39W 


7.73 


0.80 


56M 


llW 


9.89 


8.19 


10.2 


9.2 


N7537 


Sbc 


13.57W 


32.74W 


17.40W 


7.35 






0.80 


42 ra 


5W 




10.21 


10.2 


9.0 


Totals: 


86 


86 


86 


57 


76 


15 


15 


68 


61 


77 


48 


80 


80 


66 



Notes: Column (1): Galaxy ID. Column (2): Morphological type. Column (3): Galaxy B-band magnitude. Column (4): Distance 
modulus. Column (5): Nuclear star cluster magnitude. Column (6): Nuclear star cluster stellar mass. Column (7): Nuclear stellar 
disk magnitude. Column (8): Nuclear stellar disk stellar mass. Column (9) CMO mass-to-light ratio, in the band appropriate to the 
magnitudes given in Column (5) and Column (7). Column (10) Galaxy velocity dispersion. Column (11): Galaxy effective radius. Column 
(12): Galaxy dynamical mass. Column (13): Galaxy K-hand magnitude. Column (14): Galaxy stellar mass. Column (15): Spheroid 
stellar mass. 

References: fal IBinggeli et all I|1985I 1. [b] R C3. fcl ICardone fc Serenol I|2005I 1 [d] IMei et all H2007I ). [e] Assigned the median distance for 
Virgo from IMei et al l H2007I ). [f] GS09. fgl IGraharrJ J2012bl). [h] BGP07. All ma gnitu des are if -band, [i] g- band magnitudes from 
ICote" et alj l|2006l l. |i| F06. [k] IGraham et al.l ||201H ) [lj IFalcon-Barroso et at] (2003 1. [m] IFerrarese eTaTl (|2006al V 

obtained total JT-band magnitudes, m.K, for 80/86 galax- 
ies from the 2 Mic ron All Sky Survey (2MASS) Extended 
Source Catalogue (jJarrett et al.l feOOO) . To con vert to ab- 
solute magnitudes we used the distances from IMei et al.l 
(poHl for the F06 sample, and from BGP07 and GS09 
for the corresponding galaxies. We adopt an uncertainty 
of 0.25 mag for all absolute magnitudes. 

We derive dynamical masses using the simple but pop- 
ular virial estimator: Md yn — ota^H e /G, where R e is the 
effective half-light radius and a e the luminosity-weighted 
velocity dispersion measured within a 1 R e aperture. Fol- 
lowing F06, we used a value of a — 5 for the F06 galaxy 
sample. The virial factor a can take on a range of values 
(jBertin et al. 2002) depending on the radial mass distri- 
bution. By comparing virial estimator derived masses 
to the results of more sophisticated dynamical models, 
ICappellari et al.l J2006) found that, in practical situa- 
tions when working with real data, a = 5 provides a vir- 
tually unbiased estimate of a galaxy's dynamical mass 
within 1 R e .. They found this to be true for galaxies 
with a broad range of Sersic indices, n = 2 — 10 and 
bulge-to-total ratios, B/T ~ 0.2 — 1.0. They caution, 
however, that their result "strictly applies to virial mea- 
surements derived... using 'classic' determination of R e 
and L via R 1 / 4 growth curves, and with a e measured in a 
large aperture." Based on their findings we conclude that 



represent a disparate set of meas urements corrected to 
the sa me 'standard' aperture as in IFalcon-Ba rroso et al.l 
(2002]). We adopt an uncertainty of 10 % for all u mea- 
surements. 

We considered correcting the standard aperture mea- 
surements to Re measu rements based on equation 1 of 
Cappellari et al. (2006), from which we derived a mean 
correction to the a measurements reported here of 2.6 
% (ranging from +4 % to —6 % for individual objects). 
However the error on the derived correction for individ- 
ual objects is significant, ~ 10%, and much larger than 
the typical correction of ~ 2.5% for a given measure- 
ment. The correction is not correlated with host galaxy 
<T, hence will not introduce a systematic error in our un- 
corrected a values. Given this uncertainty, and that we 
were unable to derive an aperture correction for all our 
objects due to missing R e measurements, we opted not 
to apply any aperture correction. 

We obtained total apparent 5-band magnitudes for all 
galaxies following F06. For the BGP07 and GS09 galax- 
ies we obtained tub from Ide Vaucouleurs et alj (|199lL 
RC3). For the F0 6 galaxies we obtained tub from 
IBinggeli et al.l (|1985h . reduced to the RC3 system using 
the relation given in the HyperLeda database. We note 
that this ap proach fails to fully corr ect for dust in disc 
galaxies (see IGraham fc Wor lcv 2008) . We therefore also 
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the virial estimator is a reasonable approximation of the 
dynamical mass for galaxies of types Sa and earlier - we 
do not determine Mcai.dyn for galaxies of morphologi- 
cal type Sb or later, as these are heavily disk-dominated 
systems for which the virial estimator has not been cali- 
brated. 

For the BGP07 and GS09 objects we use R c values 
from the RC3 (which are determined from i? 1 / 4 curve-of- 
growth fits to the surface bright ness profile) . For the F06 
objects we use R c values from iFerrarese et al.l (|2006bT > 
which are derived from Sersic R 1 /™ fits to the observed 
surface brightness profile. For these 51 objects F06 re- 
port a range in n from 0.8 to 4.6 (with 78 per cent of 
galaxies with n in the range 1 to 2.5). We compared 
Sersic-based R e ,s for the subset of FOG galaxies for which 
RC3 RV 4 -based R c dov were available (22 objects) and 
found a one-to-one correlation, with the F06 R e ,s be- 
ing systematically 30 ± 7% larger than the R ej deV val- 
ues. For these comparison galaxies n ranged from 1.1 to 
4.6 (with 72 per cent of galaxies with n in the range 
1 to 2.5), representative of the full 51 objects. This 
suggests that, after we apply this correction to the F06 
R c ,s (Rc.doV = 0.77R 0j fo6), the use of Sersic fit based 
R e ,s will not significantly bias Mcai.dyn for the F06 ob- 
jects. While we find good agreement for this small sam- 
ple of galaxies for the specific methods used to determine 
R e by the respective authors, we caution that in gen- 
eral Sersic and R 1 / 4 based R e t ypically show significant 
differences (|Truiillo et al.l I2001D . After excluding disk- 
dominated galaxies of type Sb or later we were able to 
derive Mcai.dyn for 48/86 galaxies; all objects for which 
a a and R e measurement was available. Typical errors 
on Mcai.dyn are ~ 50 per cent. 

We additionally determine stellar masses for the full 
galaxy, Meal.*, and for the spheroidal component, 
Msph,*- To determine Meal,* we multiplied the to- 
tal galaxy luminosity of each object by the appropriate 
mass-to-light ratio. For all objects we used Mk from 
2MASS (excluding galaxies with no 2MASS M K total 
magnitude) and, following BGP07, assu med a standard 
mass -to-light ratio, M/Lk — 0.8 (see also Boll fc de Jond 
2001). Mcai,* was determined for 80/86 galaxies. All 
magnitudes an d colors were correcte d for Galactic extinc- 
tion following [Schlegeletall ([l998|). The data was not 
corrected for internal extinction, though we note that for 
most galaxies Mk and Mj Lk are minimally affected by 
dust extinction. 

Mg p h.* was determined by multiplying the total 
spheroid magnitude of each object by the appropriate 
mass-to-light ratio as described above. We use the 
spheroid masses provided by GS09 for their galaxies. We 
adopt the same spheroid masses as BGP07, obtained 
by multiplying their iiT-band spheroid magnitudes by 
M/Lft- = 0.8. For the F06 galaxies we use the galaxy stel- 
lar masses described above, excluding 16 galaxies classi- 
fied as SO or dSO as they are likely to contain a large 
scale disk and no bulge-disk decomposition is available. 
This resulted in Ms p h,* for 66/86 galaxies - all galaxies 
for which a spheroid mass or spheroid magnitude and an 
optical color were available. Typical errors on Meal.* are 
~ 25 per cent and on Ms p h,* ~ 40 per cent due to the 
increased uncertainty in separating the spheroid compo- 
nent of the galaxy's light. 



2.3. Supermassive black hole galaxy sample 

To compare to our nuclear star clust er sample, we 
take t he supermassive black hole sample of lGraham et al.l 
(|2011[ ). This sample consists of 64 galaxies with di- 
rectly measured supermassive black hole masses. The 
host gala xy velocity disper sions for this sample are pre- 
sented in IGrahamI (l2012bl). the host galaxy B- an d K- 
band luminosities in IGraham fc Scottl ([submitted!) and 
the distance to each object in IGraham et al.l (|2011l ). 

We also determine derived quantities, Mcai.dyn, Mcai,* 
and Msph,* for the supermassive black hole host galax- 
ies following the approach for the nuclear star cluster 
host galaxies described above. Briefly, we derive MGai,dyn 
from the Viria l estima tor, using the velocity dispersions 
from IGrahamI ()2012bl ) and R e from the RC3. This al- 
lowed us to derive Md yn f° r 40/64 galaxies. We de- 
rive galaxy stellar masses, Meal,* for the supermassive 
black hole galaxies as for the nuclear star cluster galax- 
ies, using the galaxy if-band magnitude and an assumed 
M/L= 0.8. We derive M Ga i,* for 59/64 galaxies - all 
objects with an available if -band magnitude. To derive 
spheroid stellar masses we make use of the spheroid mag- 
nitudes, msph presented in lMarconi fc HuritT (120031) (with 
the e xcept ion of NGC 2778 and NGC 4564, see IGrahamI 
120071) and lHaring fcRixl ([20040 These were then mul- 
tiplied by an appro priate M/L deter mined using the re- 
lations presented in lBell et alJ ((2003), with optical colors 
obtained from the HyperLeda database. We determined 
Msph,* for 39/64 SMBH galaxies - all objects with an 
available ms p h and optical color. 

3. ANALYSIS 

We use the B C ES linear fitting routine of 
lAkritas fc Bershadvl (jl996l ). which minimizes the 
residuals from a linear fit taking into account mea- 
surement errors in both the X and Y directions. We 
adopt an orthogonal minimization, BCES(Orth), 
which minimizes the residuals orthogonal to the linear 
fit. An orthogonal regression provides a symmetrical 
treatment of the data, that is, swapping x and y 
data with each other still produces the same linear 
fit. This is preferred when one is after the underlying 
physical relation, r eferred to as the "theorist's question" 
(|Novak et all 12006). To compare the two sets of scaling 
relations it is important to use the same minimization 
technique because "the different regression methods 
give different slopes even at the population level" 
(| Akritas fc Bershadvl fl996l ) . We found that, in Monte 
Carlo simulations of a mock sample of NCs and SMBHs 
drawn from a single common CMO scaling relation, 
the BCES(Orth) minimization proved most robust at 
recovering the same relation for both sets of datapoints. 
We therefore conclude that, to assess whether the two 
observed datasets are drawn from a single common 
scaling relation, the BCES(Orth) minimization is the 

3 Recent work s by IBeifiori et al.l pOLl ). ISani et al.l lf20l"lT ) and 
IVika e t al. (2012) presented new bulge-to-disk decompositions for a 
significant number of galaxies in our SMBH sample. We elect not to 
make use of these new values for now because the agreement on the 
bulge-to-total flux ratios between the three authors is poor and it 
is unclear which provides the more accurate spheroid luminosities. 
For example, the three authors find bulge- to- total ratios of 0.78, 
0.51 and 0.36 respectively for NGC 4596, with similar significant 
variations in spheroid luminosities. 
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TABLE 2 

Nuclear cluster and black hole scaling relations 





Relation 


a 


err(a) 


b 


err(b) 


C" rms 


r 


N 


Notes 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




M R + 16. 9, log M«rp 


6.44 


0.06 


-0.32 


0.05 


0.55 


-0.64 


76 






M B + 19. 9, log Mrb 


8.20 


0.09 


-0.65 


0.11 


0.39 


-0.41 


25 




Hi crn rp> 1 1 1 
jr. i^ui 1 1 1 


log (j/54.0,log Mjvc 


6.63 


0.09 


2.11 


0.31 


0.55 


0.62 


51 






log cr/224.0,log M SH 


8.46 


0.06 


6.10 


0.44 


0.47 


0.88 


64 






log M GaMy „/10 y - b ,log Matc 


6.65 


0.10 


0.55 


0.15 


0.50 


0.53 


41 


Ex. Sb and later 




log M GaM „ n /10 n - 3 ,log M SH 


8.47 


0.07 


1.37 


0.23 


0.46 


0.76 


40 


Ex. Sb and later 




+ 20.4,log Mjvc 


6.63 


0.07 


-0.24 


0.04 


0.52 


-0.69 


57 


E and dE only 




M K + 23.4,log M SH 


8.04 


0.14 


-0.48 


0.09 


0.40 


-0.70 


25 


E and dE only 


Figure [2] 


log Mc^/lO^Jog Mjvc 


6.73 


0.06 


0.80 


0.10 


0.53 


0.72 


71 




log M Ga; .»/10 n - 3 ,log M sw 


9.40 


0.32 


2.72 


0.69 


1.03 


0.55 


59 






log M Sph , t /10 a - b ,log Mjvc 


7.02 


0.10 


0.88 


0.19 


0.63 


0.64 


57 






log M Sl ,h,»/10 11 - 3 ,log M BH 


8.80 


0.11 


1.20 


0.19 


0.63 


0.65 


39 





Column (1): X and Y parameters of the linear regression. Columns (2)-(5): Slope b and zeropoint a, and their associated error, from the 
best-fitting linear relation. Column (6): Root mean square (rms) scatter in the log Mcmo direction. Column(7): Spearman r coefficient. 

Column (8): Number of data points contributing to the fit. 
Fits of the form log y = a + blogx (or log y = a + bx for Mb and Mj(), were performed using the BCES(Orth) regression. 



appropriate choice. We note that the BCES(Orth) 
method is a common technique used in determining 
linear scaling relations and has been found to pro- 
duce results co nsistent with other symmetric linear 
regressions (e.g. lHaring fc Rixl 12004). We additionally- 
calculated the Spearman's rank correlation coefficient 
for each of the correlations reported in Table [5J With 
the exception of the Mbh-M^ relation, the probability 
that the given values of r could arise if the quantities 
were not correlated is less than 0.01%. For the Mbh~Mb 
relation this probability is 0.11%. 

3.1. Nuclear star cluster mass scaling relations 

We have derived scaling relations connecting the nu- 
clear star cluster mass to various properties of their host 
galaxies: B- and JT-band luminosity, Mb and Mk, ve- 
locity dispersion er, galaxy dynamical mass MQaLdyn and 
galaxy and spheroid stellar masses, Meal,* and Ms P h.*- 
These linear scaling relations are presented in Table [21 
At present, given the relatively small samples and the sig- 
nificant errors on the derived nuclear star cluster masses, 
there is no compelling reason to fit more complicated bro- 
ken or non-linear scaling relations to our data, though 
future studies with improved data may reveal additional 
complexity. Our Figure Q] builds on Figure 2 from F06 
by presenting linear fits of Mnc against: host galaxy 
5-band magnitude M^, velocity dispersion a, and virial 
mass MGai,dyn- In Figure[2]we present fits of Mnc against 
M K , M Ga i,* and M Sp h,*. 

We find a slope of 2.11 ± 0.31 for the M NC — a rela- 
tion. This is significantly shallower than that reported 
by F06 (4.27 ± 0.61), though in better agreement with 
recently reported slopes of 1.57 ± 0.24 and 2.73 ± 0.29 
(|Grahaidl20T2at ILeigh et al.ll2012D . The principal differ- 
ence between the F06 study and the recent findings of 
a shallower slope of ~ 2 is the inclusion of nuclear star 
clusters in more massive galaxies with a > 200 km s _1 
(the F06 sample was limited to nuclear star clusters in 
host galaxies with a < 150 km s _1 ). The exclusion 
of NDs from our fits (open blue symbols in all figures) , 
which are typically an order of magnitude more massive 
than NCs, also contributes to our flatter slope, and ac- 
counts for the diff erence between our slope and that of 
ILeigh et al.l (|2012l ). 



We find a good correlation between Mnc and host 
galaxy luminosity in both the B- and if-band, with 



M 



NC °C 



L 



-0.6±0.1 

K 



We find a strong correlation of 
Mnc with Meal,* and a somewhat weaker correlation 
with Ms P h.* ~ this is consistent with the findings of 
lErwin fc Gadottil (|2012D , though we find a smaller differ- 
ence between the strength of the correlations than they 
report(0.72 and 0.65, compared to their 0.76 and 0.38). 

We find a shallow slope of ~ 0.5 for the Mnc - Moai.dyn 
relation, which is significantly flatter than the slope 
1.32 ± 0.25 reported by F06. We again attribute this dif- 
ference to the inclusion of many more massive galaxies in 
our sample (though we caution that the Virial-estimator 
based dynamical masses used here and in F06 have sig- 
nificant errors). Bearing in mind that the relations were 
not constructed to minimise the scatter in the Mcmo di- 
rection, the Mnc — MQai.dyn relation has the lowest rms 
scatter (in the vertical Mcmo direction) of any of the 
NC scaling relations (though it is not significantly tighter 
than either the Mnc — M#- or Mnc — Mc a i,* relations). 

3.2. Supermassive back hole mass scaling relations 

In this subsection we derive a set of six scaling relations 
involving black hole masses. Except for the Mbh — cr 
relation, which invol ves galaxies of all ty pes and is es- 
sentially a copy from lGraham et ail (|2011f) . due to avail- 
able data these relations predominantly involve massive 
galaxies and spheroids. As such we have not included 
the developments which reveal a bent nature to the 
other five relations at lower masses (e.g. [Graham 2012a; 
iGraham fc Scottllsubmit tcdV However, these 'bends' are 
such that the lower mass systems define steeper relations 
than shown here, which only emphasises the differences 
with the NC scaling relations discussed in the following 
section. 

We emphasise here that most of the relations we de- 
rive for supermassive black holes are for comparison only 
and do not represent the state-of-the-art in supermassive 
black hole scaling relations. For this reason we derive 
only simple linear fits to the available supermassive black 
hole sample for each variable. We do not distinguish be- 
tween barred and unbarred galaxies (see iGrahaml 120081 
for a discussion of the offset nature of barred galaxies in 
the SMBH-ct relation) , nor do we fit broken relations that 
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Fig. 1. — Mjmc and Mbh mass vs. galaxy magnitude Mg (left panel), velocity dispersion a (middle panel) and dynamical mass Meal dyn 
(right panel). Black dots indicate SMBHs, red symbols indicate NCs and open blue symbols show those objects identified as NDs. For 
the NCs and NDs the symbol indicates the sample each datapoint was drawn from: circles for F06, diamonds for BGP07 and triangles for 
GS09. The thick black and red lines indicate the best-fitting linear relations for the SMBH sample and the NC sample respectively. The 
thin dashed lines indicate the corresponding best-fitting relations from F06. A representative error bar is shown in the upper left corner of 
each pane l. We note that lower-lumino sity galaxies typically lie below the Mbh-Mb relation - consistent with the bent Mbh-Mb relation 
shown by Graham & Scott (submitted) using the spheroid luminosity. 



better describe the scaling of supermassive black hole 
mass wit h host luminosi ty (jGraham fc Scottilsubmittedl ) 
or mass (|Grahamll2012af ). Whether theMsH-Mcai,* and 
MeH-Msph,* relations are also bent is beyond the scope 
of this work but will be addressed in a future paper in 
this series. 

The linear relations we derive are presented in Table. 
[2] and are shown as the thick black lines in Figures [1] 
and [21 While we emphasise again that these linear su- 
permassive black hole scaling relations are for compari- 
son only, we briefly discuss their consistency with similar 
scaling relations presented in the literature. We note that 
our Mbh — cr relation has a slope ~ 6, whereas a slope 
~ 4 — 5 had typically been reported in the lit erature 
(|Merritt fc Ferraresd[200lt iTremaine et alJl2002t ). How- 
ever, our supermassive black hole sample now i ncludes a 
signifi cant number of barred galaxies which, as IGraharnl 
(2001) firs t observed, are offset from the unbarred re- 
lation. As Graham et al.l (|2011f) noted, including barred 
galaxies i n ones sample in creases the slope of the Mbh — & 
relation - iGraham et al.l (|2Qlll ) and IGraharnl ()2012bD re- 
port a slope of 5.95 ± 0.44 and 5.76 ± 1.54 respectively 
for their full samples of both barred and unbarred galax- 
ies, consistent with our value. While the value we re- 



port is biased by the inclusion of the barred galaxies, 
their inclusion is appropriate as we do not distinguish be- 
tween barred and unbarred galaxies in the corresponding 
Mnc — cr relation. 

The M B h - M Ga i,dyn and M B h - M Sp h,* relations 
we present, with slopes ~ 1, are typical of those re- 
ported in the lite rature (e.g. iMarconi fc Huntl [2003; 
lHaring fc~ R ix 2004])- The correlation with spheroid stel- 
lar mass is significantly tighter than with total stel- 
lar mass, consistent with othe r studies (e.g. iKormendvl 
120011 : iMcLure fc Dunlop||2002f ). Our relations are dom- 
inated by SMBHs with M B h > 1O 8 M , and IGraharnl 
(2012bJ) has shown that above this rough threshold the 
MBH/MQai,dyn ratio is fairly constant, while at lower 
masses the MBH/MGai,dyn ratio is not constant but in- 
creasingly smaller, a result which can also be seen in our 
Figure [l] where the data points fall below the extrap- 
olation of the solid black line at lower masses, causing 
the steepening of our MBH-Mcai.dyn relation. We note 
that the correlation between Mbh and host galaxy lu- 
minosity for the full sample is poor, with Spearmann r 
coefficients ~ —0.4. However, if we include only purely 
spheroidal systems (E and dE) this correlation improves 
markedly (Spearmann r = -0.81 and -0.70 in the B- and 
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Fig. 2. — Mfjc and Mbh vs - g a l ax y if-band magnitude, Mj( (for galaxies classified E and dE only, left panel), galaxy stellar mass, 
^Gal,* (middle panel) and spheroid stellar mass, Mgph,* (right panel). Colors and symbols as in Figure[T] A representative error bar is 
shown in the upper left corner of each panel. 



7<f-bands respectively), which is again unsurprising given 
that Mbh is known to correlate with the properties of 
the spheroid. 

4. DISCUSSION 
4.1. Comparison of derived relations 

In contrast to F06 and WH06 we find that NCs and 
SMBHs do not follow common scaling relations. In all six 
diagrams that we have considered, the NCs and SMBHs 
appear to follow different relations (see Table HJ) . Our 
most significant finding is that the Mnc _ c relation (with 
a slope ~ 2) is significantly flatter than the Mbh _ c re- 
lation (with a slope ~ 6 for all ga laxies, or ~ 5 for bar- 
less galaxies: [Graham et al.ll201.il ). This is in agreement 
with (Graham! (|2012bD . but in contrast to F06 who find 
an Mnc ~ o~ relation parallel to the Mbh — cr relation. 
The difference is due to: i) the exclusion of NDs from 
our NC sample; and ii) the inclusion of NCs that have 
masses higher than the SMBH/NC threshold of 10 7 M 
suggested by WH06. NDs are significantly more massive 
than NCs in comparable host gala xies and follow signif - 
icantly different scaling relatio ns ([Baleens et alJ 12007). 
iScorza fc van den Boschl (|1998[ ) showed that NDs follow 
galaxy-scale stellar disk scaling relations, extending those 
relations to much lower mass. 

In the middle panel of Figure [TJ at a CMO mass of 
~ 10 7 M , NCs are, on average, found in galaxies of 



significantl y lower a and M fl a i .d vn than SMBHs of the 
same mass. Graham (|2012a[) has revealed that the Mbh - 
Mcai.dyn relation steepens from a slope of ~1 at the high- 
mass end (Mbh ^ 2 x 10 8 Af Q ) to a slope of ~2 at lower 
masses (our slope of 1.37-1.55 is intermediate to these 
values because we fit a single linear relation to the high 
and low mass ends of what is a bent relation) . The steep- 
ening in slope at the low-mass end is in the opposite sense 
to that observed for the NCs, which have a flatter slope 
of 0.55 ±0.15. 

When considering the scaling of Mcmo with the stel- 
lar mass content of its host we find that Mnc appears 
to be driven by the total stellar mass, whereas Mbh is 
more closely associated with only the spheroidal compo- 
nent. Combining this finding with the result that Mnc 
follows much flatter relations with a and Mcai.dyn than 
Mbh does, suggests that the physical processes that lead 
to the build-up of a nuclear stellar cluster may be sig- 
nificantly different to those that drive the formation of 
supermassive black holes. A complementary view is pro- 
vided by Figure where we plot the NC mass frac- 
tion as a function of host galaxy dynamical mass and 
spheroid stellar mass, i.e. MNc/MGai,dyn vs - Mc a l,dyn 
and M NC /Ms p h,* vs. M Sp h,*. The ratio M NC /M Ga i,dyn 
shows a clear trend with Mo a i dyn, in the sense that the 
NC mass fraction decreases smoothly with MQ a ] )( jy n . We 
note that MNc/Ms P h.* is also not constant, spanning ~2 
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Fig. 3. — Left panel: Nuclear star cluster mass, Mjmc as a fraction of host galaxy dynamical mass, Mgal dyn vs. Mqj ,dyn- Right panel: 
Same as left panel except using host spheroid stellar mass, Ms p h,*- Colours and symbols as in Figure[T] The left panel shows a clear trend 
of decreasing CMO mass fraction with MQ a i i( j yn . While any trend is less clear with Mgph *, we find that the CMO mass fraction still spans 
a large range of ~ 2 orders of magnitude. In both panels nuclear stellar disks (open symbols) are significantly offset to higher CMO mass 
ratios. 



orders of magnitude from 0.02 per cent to 2 per cent. 

4.2. Galaxies hosting supermassive black holes and 
nuclear star clusters 

The GS09 sample of galaxies contain both an NC and 
an SMBH. It is likely that the SMBH and NC in these 
galaxies interacted in some way during their formation, 
hence additional physical processes may have influenced 
their scaling with their host galaxy. It is unclear what 
exactly the result of any interaction may be: it has been 
suggested that (i) the presence of a single SMBH may 
evapo rate the NC (|Ebisuzaki et alj|2001t IQ'Learv et al.l 
2006), (ii) a binary SM BH may heat and erode the NC 
(jBekki fe Grahamlfeoitil) . and that (hi) s ome SMBHs are 
in part, built up by the colli sion of NCs (|Kochanek et al.l 
U9871 iMerritt fc Poon| [2"004). Additionally, given the ex- 
istence of some dual-CMO galaxies, it is likely that some 
of the supposed NC- or SMBH-only galaxies in our sam- 
ple contain an undetected SMBH or NC, respectively. 
Because of this it is unclear whether it is correct to in- 
clude or exclude the GS09 galaxies from our main NC 
sample. More importantly, while including the GS09 
galaxies docs affect the NC scaling relations our conclu- 
sions do not depend on whether we include them. 

5. CONCLUSIONS 

We have revised three NC scaling relations and addi- 
tionally presented three new scaling relations involving 
NC mass and host galaxy properties. We have also con- 
ducted a comparison of the scaling relations for NCs and 



SMBHs for the largest sample of objects to date. Our 
principal conclusions are: 

i) The Mnc-c relation is not parallel to the Mbh - ct re- 
lation when nuclear disks are properly identified and 
excluded and recent identifications of NCs i n massive 
galaxie s are included, in agreement with iGrahaml 
(1201 2bl ). 

ii) Nuclear star clusters and black holes do not follow a 
common scaling relation with respect to host galaxy 
mass, in agreement with BGP07. 

hi) The nuclear cluster scaling relations are consider- 
ably shallower than the corresponding supermassive 
black hole scaling relations. This is true for the re- 
lations involving host galaxy: a, luminosity, dynam- 
ical mass and stellar mass. 

iv) The dominant physical processes responsible for the 
development of NCs and SMBHs, in relation to their 
host galaxy or spheroid are suspected to be different 
given the above findings. 

v) The NC mass fraction, with respect to the mass of 
its host galaxy or spheroid, is not constant, span- 
ning ~2 orders of magnitude. The NC mass fraction 
decreases in more massive galaxies. 
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